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O \ Abstract 

We show that the infinite-dimensional representation of the re- 
cently introduced Logistic algebra can be interpreted as a non-trivial 
qq . generalization of the Heisenberg or oscillator algebra. This allow us to 

C*~) ! construct a quantum Hamiltonian having the energy spectrum given 

' by the logistic map. We analyze the Hamiltonian of a solid whose 

■ collective modes of vibration are described by this generalized oscil- 

"O lator and compute the thermodynamic properties of the model in the 

two-cycle and r = 3.6785 chaotic region of the logistic map. 
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1 Introduction 



In the last years, complex systems have attracted a lot of attention. In par- 
ticular, there is an intrinsic theoretical interest in constructing a Hamiltonian 
system having an energy spectrum that is or quasi-periodic or self-similar, 
and/or chaotic Enhancing the interest in describing such Hamiltonian 
system is the fact that some models on quasicrystals have a quasi-periodic or 
fractal energy spectrum p], || |, §, |], [/], [8j . On the other hand, one paradig- 
matic example of map that exhibits some of these features is the logistic map. 
As is well-known, this map describes at the Feigenbaum point an example of 
a fractal system and appearing after this point a chaotic region with chaotic 
bands and self-similar patterns [Q. 

Recently, it was developed a three-generator algebra, called Logistic al- 
gebra ||10|| , where the eigenvalue of one generator is given by the logistic map 



. We show that the infinite-dimensional representation of this algebra can 
be interpreted as a non-trivial generalization of the Heisenberg or oscillator 
algebra and call the associated oscillators, logistic oscillators. 

We use these logistic oscillators to construct a quantum Hamiltonian, 
which is a generalization of the quantum harmonic oscillator, that has the 
energy spectrum described by the logistic map. We apply these ideas to 
construct a Hamiltonian describing quasi-particle vibrations of a solid with 
N atoms where each quasi-particle oscillates as a logistic oscillator. 

In section 2, we discuss the logistic algebra and its interpretation as gen- 
eralized Heisenberg algebra, in section 3 we construct a model for a solid 
where the collective modes of motion are described as logistic oscillators and 
compute the thermodynamic functions of the model in the two-cycle and 
r = 3.6785 chaotic region of the logistic map x n+1 = rx n {l — x n ) . Section 
IV is devoted to our conclusions. 



2 Algebraic origin of the model 

The model we are going to discuss in sections 3 and 4 has its origin in an 



algebraic structure called Logistic algebra . In this section we present the 
Logistic algebra and show that this algebra can be interpreted as a non-trivial 
extension of Heisenberg algebra. 

Let us consider the algebra generated by J , J± , described by the rela- 
tions P| 
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JiJ + = J + J i+1 z = 0,l,2,... , (1) 

J-Ji = Ji+iJ- , (2) 
J + J- — J-J + = —a(J — Ji) , (3) 

where J_ = J+ , j/ = Jj and "a" is a real constant. Moreover, 

J m = rJi(l-Ji) % = 0,1,2,... , (4) 

with < r < 4 . 

The hermitian operator J can be diagonalized. Consider the state |0 > 
with the lowest [] eigenvalue of J 

J |0 >= a |0 > . (5) 

Note that, for each value of we have a different vacuum that for simplicity 
all of them will be denoted by |0 >. We choose < a < 1 because with 
this condition all the future iteration will remain in this interval and the 
connection with the chaotic concepts is straightforward. Also, the allowed 
values of «o depend on r and a. Since, by hypothesis, «o is the lowest J 
eigenvalue we must have, 



J_|0>=0 . (6) 

Following the usual steps for constructing (now from lower to higher eigen- 

using the algebraic relations ex- 
|||)]) we obtain: 



values) SU(2) algebra representations [|Tl 
hibited in eqs. an d taking into account eqs. 



Jo\m > = a m \m > , (7) 
J + \m> = N m \m + l> , (8) 
J_|m + 1> = N m \m > , (9) 

where 

1 Duc to the use of the logistic map, depending on the values of r and ao considered, 
|0 > can be the state with highest weight. We emphasize in this paper the case where 
|0 > is a lowest weight vector since it is the situation considered in the following sections. 

2 Note that if we put m = —1 in the eq. (|J) we obtain consistently eq. (0). 
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N m = \Ja(a - oim+i) , (10) 

and 

a m+1 = ra m (l - a m ) . (11) 

Note that the states \m > , m > 1 , are defined by the application of J + on 
\m — 1 > . Moreover, from eqs. (|7| - 1^) we easily obtain a general expression 
for \m >, 

\m>= f ^ TW (J + r\0> ■ (12) 

Of course, since the eigenvalues of Jo are given by the logistic map (eq. 
flll|)), their values as m increases can have an irregular behavior depending 
on the values of r and «o , and the dimension of the representation. Note 
that, unlike 577(2) algebra where the states obtained by the application of 
J + always have higher J eigenvalues, for the logistic algebra this depends on 
what values of r and «o we consider and the level of iterations (the number m 
of \m >) we are. For instance, for r = 3 and a m = 0.5 we have a m+ i = 0.75 
, i.e., J + rises the Jo eigenvalue of \m > . On the other hand, for r = 1.5 and 
a m = 0.5 we have a m+ i = 0.375 and in this case J + lowers the J eigenvalue 
of \m > . Moreover, due to the non-regular behavior of the logistic map, it 
may happen for J + that even having started as lowering the Jo eigenvalue of 
\m > it rises the Jo eigenvalue of J + |m > for a given level m of iteration of the 
logistic map. For instance, for r = 2.75 and a m = 0.9 we have a m+ i = 0.247 
and a m+ 2 = 0.5122 . 

Let us now consider the operator 

C = J + J_ + aJ = J- J+ + aJi . (13) 
Using the algebraic relations (eqs. ([l],Q,[3D) it is easy to see that 

[C,J } = [C,J ± ] = , (14) 
i.e., C is the Casimir operator of the algebra. In fact, we arrive easily at 

C\m >= c \m > , (15) 
with Co = aao independent of m . 
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With respect to matrix representations of the Logistic algebra there are 
finite-dimensional matrix representations corresponding to the n-cycle solu- 
tions of the logistic map and infinite-dimensional ones relative to the n-cycle 
and to the chaotic regime of the logistic map. Here we present some exam- 
ples: 

1. Two-Dimensional Representations 



Jo 



Q 





J4 





N 



J_ = Jl 



The allowed values of r and ao are determined by the equation N± 
that ]Vq 7^ . There are two non-trivial solutions, 



(16) 



such 



2r - 3 



2r 



The solution gives > (%i implying a > 



, while a 



(17) 

< cti gives a < 0. 



For both cases r > 3. We will use this solution in the next section. 
2. Three-Dimensional Representations 



Jo 



f a 

«! 

V a 2 
















J + = 




No 
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Ni 






The allowed values of r and «o are computed from N 2 
3. Infinite-Dimensional Representations 



j_ = 4 • (18) 
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< Oto 
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V 



J_ = Jl 



J 



(19) 

The allowed values of r and «o can be computed for instance for a < from 
Nm = H(a< m+ i — «o) by imposing a m > «o for all values of m > 1 . In fig. ^ 
we show a half-leaf region with the allowed values of r and «o satisfying the 
above requirements. These solutions will be used in the following section. 
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Figure 1: Region of allowed values for ao an d r - 



Let us show now an interesting connection of this algebra with the Heisen- 
berg algebra. The Heisenberg algebra is generated by the elements A and 
A + satisfying the relations 

AAt-AU = 1 , (20) 
NA^-A^N = A^ , (21) 

with iV = A^A is the number operator. Note that eqs. fll], @), for i = can 
be seen as defining equations for J\. The Heisenberg algebra comes naturally 
if we put in eqs. (|l]fj, |3|) J_ = A, J + = A', Jo = N, J\ = Jo + 1 and a = — 1 
. It can be easily verified that we do not have in this case finite dimensional 
representations and the Casimir operator is identically null. 

In summary, Heisenberg algebra is the special case of the defining relations 
given by eqs. (|l],|2|, [|) where instead of taking the relation given by eq. (01) we 
consider the simpler one J\ — Jo + 1 . In other words, the Logistic algebra 
can be interpreted as an extension of the Heisenberg algebra where instead 
of the simple iteration J, + i = J, + 1 we take the logistic map for Jj + i as in 
eq. (U) . Clearly, it is also possible here to consider other maps; this study 
is under progress . 
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Of course, since the Heisenberg algebra is a master algebra in physics, it 
is a natural step to investigate the possible consequences of the logistic gen- 
eralization, explained before, in physical problems. In the following sections 
we apply this generalized Heisenberg algebra to a collective modes of motion 
of N atoms. 



3 Model and thermodynamic properties 

Let us consider the Hamiltonian of a quantum system of quasi-particles 
described by N independent, localized, "oscillators" of the form 

N 

H = T,e q f , (22) 

q=l 

where {Jq} is a collection of N independent oscillators, each of them de- 
scribed by the algebra (1-3), and e q is a parameter associated to the energy 
of the q-th oscillator. We are then considering independent collective exci- 
tations with a non-trivial spectrum specified by the eigenvalues of Jq. For 
the solution J_ = A , J + = , J = iV , J i+ i = Jj + 1 and a = — 1 of 



the algebra (1-3), the Hamiltonian fl22 ) describes the well-known system of 



N independent, localized, harmonic oscillators. On the other hand, by con- 
sidering the logistic generalization, eq. (|22|) becomes the Hamiltonian of a 
system of quasi-particles described by N independent, localized, logistic oscil- 
lators. We interpret Jl, J\ and Jq as annihilation, creation and generalized 
number operator, respectively, of the g-th oscillator. Note that the energy of 
the g-th oscillation- mode in a state \m > is given by the product of e q times 
the eigenvalue of {Jq} applied on that state. The eigenvalue ofi m indicates 
that the g-th oscillation-mode is in the state \m >. We are adopting this 
model due to its simplicity but these logistic oscillators could also be used in 
more complicated models as for example in disordered systems. 
The partition function of the model (|22D 



Z = Trexp(-pH) , (23) 

with (3 = (ksT) -1 and ks the Boltzmann constant, factorizes into a product 
of single particle partition functions, 

Z = I[Z q , (24) 
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Z q = £exp(-/3e 9 a£J , (25) 



m=0 



where the trace was performed using the basis described in eqs. (|| - [L4D 



and a m+ i = ra m (l — a m ) . We take the simplest case where a and r are 
independent of q. 

We suppose that the dispersion relation of the quasi-particle is given by 
(equivalent to the Debye approximation), 

e <? = e(?) = 7<? , (26) 

and we enclose the system in a large 3-dimensional volume V. Replacing, in 
the usual way (since we are considering phonons with a spectrum different 
from the harmonic oscillator one), the sum over particles by an integral over 
a g-space, 

E A, (27) 

we obtain, for the logarithm of the partition function, after integrating over 
the angular variables, 

V fiu (™ \ 

InZ = — / dqq In > exp(-[3~fqa m ) , (28) 
2vr 2 Jo Vm^o / 

where this integral is evaluated over a finite q-range corresponding to a finite 
number of oscillators, and qu is the larger possible number q. The mean 
energy of the solid becomes after defining a new variable rj = fi^yq, E = ■yqM, 
T = E /k B and A = X^f 

E = - 9 -^ = AEo&f [ T0/T dW ^r^f^ ^ • (29) 
d(3 V T ; Jo E^=o ex P(-^«m) 



Let us study the integrand of eq. ( p9|) . The sum is performed over the 
integer m that corresponds to the level of iteration of the logistic map since 
a m is given by this map. In what follows we shall consider two cases: an 
example of the two-cycle and another one corresponding to the chaotic region 
of the logistic map. 

At a given approximation, in the two-cycle region of the logistic map (3 < 
r < 3.449489...), the iteration runs over transient states before reaching 
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the asymptotic two levels, which are infinitely degenerated. Clearly, when 
the degeneracy g of the two levels goes to infinity the contribution of the 
transient states disappears and only the contribution of the states related 
to the asymptotic levels remains. The measure is concentrated on the two 
asymptotic levels. The effective expression for the energy in the infinite g 
limit is given by: 

1q Jo exp{—r]a ) + exp(— r]a + ) 

For the specific heat at constant volume we have: 



C v 
where 



'dE? 
dT 



Ak„(L? 



v 



T /T 



dw 3 f(v) 



Jh4 

T 



/(T /T) 



(31) 



a 



exp(— rja ) + a + exp(— r/a + ) 



(32) 



exp(— rja ) + exp(— r]a + ) 

In figure |2] we display e = E/AEq times t = T/To; in figure |3] we show 
C = Cy/Co times t with Co = Aks ■ These are typical graphics for two- 
level systems since after the transient states what remains is the two-cycle 
situation. For higher-cycle regions of the logistic map we shall have the 
typical behavior of a system with a finite number of levels. 

If we calculate the entropy from eq. fl2"gp we see that it diverges, since 
the degeneracy factor g goes to infinity. The renormalized entropy Sr = 
(S/k — (A/3) log g) I A can be calculated and expressed as: 



S 



r 



T 
% 

T /T 



rT /T 

J drjrf log (exp(— r/a + ) + exp(— r]a~)J + 



dr]r] 3 f(rj) 



(33) 



More interesting is the behavior of the system we are analysing for the chaotic 
region. In this case we have as before transient states with the difference that 
instead of having a finite number of asymptotic levels we have a continuum 
of levels similar to the classical continuum levels in a classical system. Thus, 
after dropping the transient states, as the measure is concentrated on the 
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Figure 2: Two cycle energy versus temperature. Continuous line, r = 3.1. Broken 
line, r = 3.35 

chaotic region, the system is better described by a density function that 
represents the number of hits of the logistic map in the interval [0, 1] . In 
this case the mean energy is given by 



T r T o/ T 
E = AEni—) 4 I dr] V 3 F(r)) 



o Jo 



(34) 



with 



F(v) 



J? dfttf (fl exp(-itf) 



(35) 



fidZH(Z)exp(-vt) ' 
where H(£) is the density function. The specific heat at constant volume 
becomes 



'dET 
df 



J-Cl 



■To /T T T 

I / q d VV 'F( V ) - (^) 4 f(^; 



(36) 



In figures f| and § we show E/AEq and Cv/Cq times T/Tq in the chaotic 
region for r = 3.6785 . 
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Figure 3: Specific heat for a two cycle. Continuous line, r = 3.1; broken line, 
r = 3.35. 

Figures ^ and [| exhibit the typical low-temperature behavior of a classical 
system as we had already anticipated, since we have a continuum energy level 
in the chaotic band. As the spectrum is limited from above, the behavior of 
the specific heat, at high temperatures, for any value of r, is proportional to 
1/T 2 as expected from systems that present the Schottky anomaly. 

We also show in fig. |6| the normalized histogram and the density function 



[ ( (tt/2)^/(£ - 0.266) (0.726 - £) ) if 0.266 < f < 0.726 , 

H ® = \ ( (vr/2)^(£ - 0.728)(0.922 - f) ) _1 if 0.728 < f < 0.922 , 
[ otherwise , 

(37) 

we used in order to compute the mean energy, the specific heat and the 
entropy. 

After dropping the transient states, the entropy for the chaotic region 
S c haos = S/Sq, where S = C = Ak B , can be expressed as: 
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Figure 4: Energy versus temperature for a chaotic spectrum. 



. 3 



0.2 



. 1 




0.5 1 t„ 

Figure 5: Specific heat versus temperature for the chaotic spectrum. 



s, 



chaos 



T 
% 

T /T 



| Q T ° /T d Vr] 2 log QT 1 d£#(£) exp(-^)) + 



drjrfFirj) 



(38) 



The entropy of the chaotic band, showed in fig. [7|, presents a curious behavior. 
In fact, its low-temperature behavior is typical of a classical system, with a 
negative divergence as T — ► . On the other hand, its high-temperature 
behavior is typical of a system with a limited spectrum, found mainly in 
quantum systems. The origin of this hybrid behavior is the fact that the 
neighbor levels in the chaotic band have no minimum distance among them, 
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.005 




Figure 6: Histogram of the logistic map for r = 3.6785 and the function we used 
(full line) to calculate the thermodynamic functions. 



since they are dense inside the chaotic band. This is somewhat equivalent 
of taking h — > limit, thus reobtaining classical low-temperature behaviors. 
Note that in this case this limit is not imposed, but it is intrinsic of the 
system since the commutation relations are always different from zero. 
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Figure 7: Entropy versus temperature for a chaotic spectrum. 



4 Conclusion 



We construct, based on an algebra developed in ||10|| , a Hamiltonian of a 
quasi-particle that presents an energy spectrum whose energy levels are gen- 
erated by the logistic map. Depending on the parameter r of the logistic map, 
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the energy levels can be finite (corresponding to cycles of the logistic map) 
or chaotic (corresponding to the chaotic bands of the map). We study the 
thermodynamic properties of a Debye-like solid constituted by these elemen- 
tary quasi-particles and we exhibit the behavior of some thermodynamical 
quantities like internal energy, specific heat and entropy These functions, as- 
sociated to the chaotic spectrum, present a mixed aspect, with both classical 
and quantum-like typical behaviors. This is consequence of the fact that the 
spectrum in the chaotic region is continuous, similar to the spectra of clas- 
sical systems, with no separation between neighboring levels. On the other 
hand, the thermodynamic quantities related to the cycles are equivalent to 
systems with a finite number of energy-levels. 

It is interesting to note that the algebraic formalism developed in section 
H works consistently for a large class of maps Jj+i = f(Ji) ■ Of course, 
changing eqs. (|J |TT| ) implies a different representation theory of the algebra 
(|l|-|3|) and a different physical Hamiltonian. 

A classification of the analytical functions / under a stability theory would 
lead us to determine the different Hamiltonians associated to the different 
kinds of attractors of the map /. A systematic study of different non-trivial 
relations Jj+i = f(Ji) and their consequences on the Hamiltonian spectra is 
under study. 

We acknowledge interesting remarks from C. Tsallis. We also thank fi- 
nancial support by CNPq and Pronex/MCT. 
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